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Abstract 

A new approximation formalism is applied to study the bound states of 
the Hellmann potential, which represents the superposition of the attractive 
Coulomb potential —a/r and the Yukawa potential 6exp(— 5r)/r of arbitrary 
strength h and screening parameter 6. Although the analytic expressions for 
the energy eigenvalues E^ i yield quite accurate results for a wide range of 
n,£ m the limit of very weak screening, the results become gradually worse 
as the strength b and the screening coefficient 6 increase. This is because 
that the expansion parameter is not sufficiently small enough to guarantee 
the convergence of the expansion series for the energy levels. 
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I. INTRODUCTION 



Adamowski [1] has presented a study of the systematics of the energy eigenvalues of the two- 
particles interacting via the superposition of the Coulomb and Yukawa potential (SCYP): 

V {r) — —a/r -\-hexjp{—5r)/r, (1) 

where a and h are the strengths of the Coulomb and the Yukawa potentials, respectively, 
and 5 is the screening parameter. It is assumed that a and 5 are positive whereas h can be 
positive or negative. The potential in (1) with h positive was first suggested by Helmann [2,3] 
many years ago and henceforth this potential will be referred to as the Hellmann potential 
irrespectively of the sign of h. The Hellmann potential has been used by various authors 
to represent the electron-core [4,5] or the electron-ion [6,7] interaction. Varshni and Shukla 
[8] used this model potential for alkali hydride molecules. Das and Chakravarty [9] have 
proposed that such a potential is suitable for the study of inner-shell ionization problems. 

The bound-state energies of the Hellmann potential for various sets of values of h and 5 
have been studied elaborately by Adamowski in a variational framework using ten variational 
parameters. The energy eigenvalues have been predicted very accurately but the calculations 
involve extensive computational time and effort. Moreover, compact analytic expressions far 
the energy eigenvalues are not obtainable. On the other hand, Dutt et al [10] have also been 
investigated the bound-state energies as well as the wave functions of this potential using 
the large- expansion technique. 

In this paper, we study the bound-state properties by applying a new methodology 
[11] based on the decompose of the radial Schrodinger equation into two pieces having an 
exactly solvable part and an additional piece leading to either a closed analytical solution or 
an approximate treatment depending on the nature of the additional perturbed potential. 
The application [12,13] of this novel treatment to the different problems in both, bound and 
continuum regions, have been proven the success of the formalism. We demonstrate here 
how such interaction potential can be simply treated within the framework of the present 
formalism. 
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One of the prime motivations of the present study is to explore the regions of vahdity of 
this approximation formahsm for the superposition of potentials such as the one in (1) which 
manifests different structures for various range of values of b and S. Our calculations reveal 
that the degree of accuracy of the predicted eigenvalues vaires appreciably for different range 
of values of b and 5 and for different quantum states. These observations have relevance in 
the context of applications of this novel treatment to realistic problems of atomic physics. 
The other motivation is that potential (1) with a — and b — —aZ can be reduced into the 
static screened Coulomb potential (SSCP) of the simple form: 

V{r) = -{aZ)exp{-Sr)/r, (2) 

where a — (137.037)"^ is the fine-structure constant and Z is the atomic number, is often 
used to compute bound-state normalizations and energy levels of neutral atoms [14,15,16,17] 
which have been studied over the past years. It is known that SSCP yields reasonable results 
only for the innermost states when Z is large. However, for the outermost and middle atomic 
states, it gives rather poor results. Although the bound- state energies for the SSCP with 
Z — 1 have been studied [18] in the hght of the shifted large- A?" method. Recently, we have 
also investigated that this novel treatment is useful in predicting bound-state energy levels 
of light to heavy neutral atoms [17]. 

The contents of this paper is as follows. In section II we breifly outline the method 
with all necessary formulae to perform the current calculations. In section III we apply 
the approach to the Schrodinger equation with Hellmann potential and present the results 
obtained analytically for the bound-state energy values upto third perturbation energy shift. 
Section IV contains the numerical results. Finally, in section V we give our conclusions. 

II. THE METHOD 

For the consideration of spherically symmetric potentials, the corresponding Schrodinger 
equation, in the bound state domain, for the radial wave function reads 



3 



2m V'n (r) 



V{r) - E^, 



(3) 



with 



V{r) 



Voir) + 



+ 1) 



2m 



+ AV{r), 



(4) 



where AV{r) is a perturbing potential and ^n(^) = Xn{r)un{r) is the full radial wave func- 
tion, in which x„(r) is the known normalized eigenfunction of the unperturbed Schrodinger 
equation whereas Un{r) is a moderating function corresponding to the perturbing potential. 
Following the prescription of Refs. [11,12,13], we may rewrite (3) in the form: 



(5) 



2m \Xn{r) Un{r) Xn{r)un{r) ^ 

The logarithmic derivatives of the unperturbed Xn{r) and perturbed Un{r) wave functions 
are given by 



^Xn{r) 



and AWn = 



(6) 



which leads to 

x:{r) 

2m Xn{r) 



W'(r) - 



n 



/2m 



--W(r) 



Voir) + 



e{i + 1) 

2m 



- Sr. 



(7) 



where £„ is the eigenvalue for the exactly solvable potential of interest, and 



2m \Un{r) Xnir)un{r)^ 



^W^ir) - -L 
V2m 



A<(r) + 2Wn{r)AWn{r) = AV{r) - Asn, 



(8) 



in which A£„ = E^^ + E^^ + E^^ + • • • is the correction term to the energy due to AV{r) 
and En = en + Ae„. If Eq. (8), which is the most significant piece of the present formalism, 
can be solved analytically as in (7), then the whole problem, in Eq. (3) is 

h 



[Wn{r) + AWn{r)Y - ^= [Wn{r) + AWn{r)]' = V{r) - E^, 



'2m 



(9) 
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which is a well known treatment within the frame of supersymmetric quantum theory 
(SSQT) [19]. Thus, if the whole spectrum and corresponding eigenfunctions of the un- 
perturbed interaction potential are known, then one can easily calculate the required super- 
potential Wn{r) for any state of interest leading to direct computation of related corrections 
to the unperturbed energy and wave function. 

For the perturbation technique, we can split the given potential (3) into two parts. The 
main part corresponds to a shape invariant potential, Eq. (7), for which the superpotential is 
known analytically and the remaining part is treated as a perturbation, Eq. (8). Therefore, 
it is obvious that Hellmann potential can be treated using this prescription. In this case, the 
zeroth-order term corresponds to the Coulomb potential while higher-order terms consitute 
the perturbation. However, the perturbation term in its present form cannot be solved 
exactly through Eq. (8). Thus, one should expand the functions related to the perturbation 
in terms of the perturbation parameter A, 

oo oo oo 

AV{r;X)=Y.\,V,{r), AW^„(r; A) = ^ W«(r), E„{X)=Y.\E^\ (10) 

1=1 i=l i=l 

where i denotes the perturbation order. Substitution of the above expansions into Eq. (8) 
and equating terms with the same power of A on both sides up to O(A^) gives 

V2m dr 

Wi'\r)Wi'\r) + 2Wr.{r)Wl^\r) - = V,{r) - (12) 

2 [W^{r)Wl^\r) + Wi'\r)Wi^\r)] - ^"^^^^"^ = Vs{r) - (13) 



2m dr 

(14) 

Hence, unlike the other perturbation theories, Eq. (8) and its expansion, Eqs. (11)-(14), 
give a flexibility for the easy calculations of the perturbative corrections to energy and 



wave functions for the nth state of interest through an appropriately chosen perturbed 
superpotential. 



III. APPLICATION TO THE HELLMANN POTENTIAL 

Considering the recent interest in various power-law potentials in the literature, we work 
through the article within the frame of low screening parameter. In this case, the Hellmann 
potential can be expanded in power series of the screening parameter 5 as [13,17,20] 

V{r) = ~ + ^- exp(-5r) = + ^ E V,{Sry, (15) 

where the perturbation coefficients Vi are given by 

V, = +1, V2 = -1/2, Vs = 1/6, = -1/24, • • • . (16) 

We now apply this approximation method to the Hellmann potential with the angular mo- 
mentum barrier 



a b , ^ , i{e+i)h 



2 



V{r) — 1 — exp(— (5r) 



Voir) + 



+ 1)^' 



2mr2 



+ AV{r), (17) 



r r 2mr2 
where the first piece is the shape invariant zeroth-order which is an exactly solvable piece 
corresponding to the unperturbed Coulomb potential with Vo(r) = —{a—b)/r while AV{r) = 

—bS + {bS'^/2)r — {bS^/Gy^ + (6(5^/24)r^ is the perturbation term. The literature is rich 

with examples of particular solutions for such power-law potentials employed in different 
fields of physics, for recent applications see Refs. [21,22]. At this stage one may wonder 
why the series expansion is truncated at a lower order. This can be understood as follows. 
It is widely appreciated that convergence is not an important or even desirable property for 
series approximations in physical problems. Specifically, a slowly convergent approximation 
which requires many terms to achieve reasonable accuracy is much less valuable than the 
divergent series which gives accurate answers in a few terms. This is clearly the case for 
the Hellmann problem [23]. However, it is worthwhile to note that the main contributions 
come from the first three terms. Thereby, the present calculations are performed up to the 
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third-order involving only these additional potential terms, which suprisingly provide highly 
accurate results for small screening parameter 5. 



A. Ground State Calculations {n = 0) 



In the hght of Eq. (7), the zeroth-order calculations leading to exact solutions can be 
carried out readily by setting the ground-state superpotential and the unperturbed exact 
energy as 



Wn=o (r) 



h £ + 1 



+ 



m {a — b) 



2 {£+l)h' 



m{a — by 



2r(n + £+l)2' 



n = 0,1,2, 



(18) 



and from the literature, the corresponding normalized Coulomb bound-state wave function 
[13,17,24] 



xPir) = exp [-/3r] x +^ [2/3r] , 

2m{a—b) 



in which N, 



and (x) = E:;=o(-1'™ 



e+i 



(19) 



1^^)]/ n\n+2e+iy. ^ normahzed constant, (3 = ^^f^ 
(n-rn)!(Tn+fc)!m.! -^'" associate Laguarre polynomial function [25]. 
For the calculation of corrections to the zeroth-order energy and wavefunction, one needs 
to consider the expressions leading to the first to the third-order perturbation given by Eqs. 
(11)-(14). Multiplication of each term in these equations by Xn(^)' bearing in mind 
the superpotentials given in Eq. (6), one can obtain the straightforward expressions for the 
first-order correction to the energy and its superpotential: 



Ei 



1) 2 

J —oo 



'b6^ 



dr, Wi'') 



'2m 1 



and for the second-order correction and its superpotential: 



dx, (20) 



.^r^-iy« (r)W^«(r) 



dr, 



/2m 1 



E(2) + w^(i) (r)l^«(x) + ^a;^ 



dx, 



(21) 
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and for the third-order correction and its superpotential: 



= r Xlir) 

J —oo 



24 



dr, 



(22) 



for any state of interest. The above expressions calculate W^^\r), H^^^^(r) and W^^\r) 
exphcitly from the energy corrections E^^ and E^^ respectively, which are in turn 

used to calculate the moderating wave function Un{r). 

Thus, using Eqs. (20)-(22), one finds the zeeroth order energy shift and their moderating 
superpotentials, for a ^ b, as 

hbNoS^ 



4(a — b)m 



2V2m{a-b) ' 



(2) ^ h%N^ jbN^ - 3L + 1) _ WiV,^ (5iVg - 3L + 1) 
° 12(a-6)2m2 16{a-b)''m^ 



hNo {a-b)mr + n'NoNi 3h%'^N^5 + 4mb{a - b) 



24V2mm2(a - by 



(3) _ h%N^ {5N^ - 3L) (5iV2 _ 3L + 1) ^^ft^^^ (5^2 _ 3/, + 1) (giya _ 5/,) 



£/0 — 



96 (a — bYrrfi 



5* + 



48 (a - 6)5m4 



^5^2 _ 3L + 1) (97V2 _ 



8\ 



(23) 



64(a — bynv' 

where A^o = (-^ + 1) , -^1 = (^ + 2) and L — £{£ + !). Therefore, the analytical expressions for 
the lowest energy and full radial wave function of the Hellmann potential are then given by 



E^=o,i = E^nlo -bS + Ej,'^ + + ^(3) ^ . . . ^ ^n=oAr) ^ xi=l^(r)wn=o/(r), 



Un=oAr) ~ exp ( f {W^'^ {x) + W^'^ {x)) dxj . 



(24) 



Hence, the explicit form of the full wave function in (24) for the ground state is 
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'0n=O,^(r) 



2m(a — b) 



{a-b)m 

exp(P(r)), 



(25) 



with P(r) = J2i=2Pi''^^ is ^ polynomial of third order having the following coefficients: 



P2 



_bNoP_ 
4(a - b) 



Nih^c 
m 



, P3 = lbc6^, 
6 



(26) 



in which 



12m(a - 6)3 



3h'bN^6 + 4m(a - by 



(27) 



B. Excited state calculations (n > 1) 

The procedures applied in the calculations of the ground states becomes extremely cum- 
bersome in the description of radial excitations when nodes of wavefunctions are taken into 
account, in particular during the higher order calculations. Although several attempts have 
been made to bypass this difficulty and improve calculations in dealing with excited states, 
(cf. e.g. [26], and the references therein) within the frame of supersymmetric quantum 
mechanics. 

Using Eqs. (6), (18) and (19), the superpotential Wn{r) which is related to the excited 
states can be readily calculated by means of Eqs. (20)- (22). Therefore, the energy shifts in 
the first excited state, with a ^ b, are: 



4m(a-6) 2V2m(a - 6) 



(2) ^ h%Nf{5Nf-3L + l) _ h%^N^{5Nf-3L + l) 
' 12(a-6)2m2 16{a - b)'^m^ 



hNi {a-b)mr + ti'NiN2 3h%'^NfS + Amb{a - b) 



24V2mm2(a - by 



5' 

— r. 



^(3) _ h%N^ {5N^ - 3L) (5iVf - 3L + 1) h%''Nt (5A^f - 3L + 1) {9Nl - 5L) 
^ ~ 96{a-b)^m^ ^ 48(a - 6)5m^ 



, h'^Nf (57Vf - 3L + 1) (QATf - 5L) 
^ 64(a - bym^ ^ ' 

Therefore, the approximated energy value of the Hellmann potential corresponding to the 
first excited state is 

En=i,e - Ei'^ -bS + e['^ + £;f ^ + ^ + • • • . (29) 

The related radial wavefunction can be expressed in an analytical form in the light of Eqs 
(20), (21) and (24), if required. The appromation used in this work would not affect con- 
siderably the sensitivity of the calculations. On the other hand, it is found analytically 
that our investigations put forward an interesting hierarchy between W^^\r) terms of dif- 
ferent quantum states in the first order with a ^ b after circumventing the nodal difficulties 
elegantly, 

»f\r) = ^ r, (30) 

2v2m(a — b) 

where Nn — {n-\-l-\-l). Therefore, for the second excited state (n = 2) leads to the first-order 
energy shift and superpotential 

' 4m(a-6) ^ ^ ' 2^2^(a - 6) ^ ' 

where N2 = (£ -|- 3) . Thus, the use of W2^^{r) in Eq.(21) gives the energy shift and super- 
symmetric potential in the second-order with a ^ b as 

(2) ^ n%Ni{5Nl - 3L + 1) _ h%'Nli5Ni-3L + l) 
' 12(a-6)2m2 16{a - b)''m^ 



{a - b)mr + n'N2N^ 2>h%'^N^5 + 4m6(a - b) 



24V2m(a - bYm? 
Finally, we obtain the third-order energy shift as 



-r. 



(32) 



n%Nl {5Ni - 3L) (5iV| - 3L + 1) h^b'^Nl (5A^| - 3L + 1) (9A^| - 5L) 



96(a — 6)^m^ 



48(a - 6)5m^ 
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h'%^Nl {5Ni -3L + 1) {ml - 5 L) 
64(a — bym^ 



(33) 



Therefore, the approximated energy eigenvalue of the Hellmann potential corresponding to 
the second excited state is 



En=2,i = 4°^ -bS + + + E)^> + 



(1) . p(2) . p(3) 



(34) 



In general, from the supersymmetry, we find out the nth state energy shifts together with 
their supersymmetric potentials with a ^ b as 



4(a -b)m ' ^ ^ ' 2x/2m(a - b) 



r. 



12(a - bfm? 



16(a - bYrn? 



a - b)mr + fi^N^^Nr^+i 3h%'^N^S + Amb{a - b) 



5^ 



2AV2rnm'^{a-by 



(3) ^ h%N'j5N', - 3L] [57V2-3L+1] W7V„^ [5jV^ - 3L + 1] [9iV„^ - 5L] 
96(a-6)3m3 ^ 48(a - 6)5m4 

^^1063^6 (5^2 _ 3^ + 1) _ 5L) ^6 



5^ 



64(a — bym^ 

where A^n+i = {n + l + 2). Consequently, the total energy for the nth state is 



(35) 



(36) 



IV. NUMERICAL RESULTS 

For the numerical work, we take a = 2 and thus our b is to be identified as the corre- 
sponding rescaled parameter in Adamowski's paper. Consequently, our energy eigenvalues 
are like those obtained by Adamowski. In Tables 1 and 2, we compute the binding energies, 
—En^i, of the lowest-lying states (from Is up to 4/) for various values of 6 = ±1, —2, —4, —10 
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as functions of the screening parameter 5 obtained from the analytic expressions given in 
Eqs. (18), (35) and (36). The dependance of the energy levels on b is shown in Table 3 for the 
states Is up to 3d. The results for the higher excited states (from 5s up to 7i) are presented 
in Table 4. The energy eigenvalues Is — 6h for the attractive Yukawa potential with S — 0.1 
as functions of b are shown in Table 5. The predicted results are then compared with the 
accurate energy eigenvalues [1] obtained by Adamowski using a high precision variational 
technique. 

Therefore, we display our results in Tables 1 and 2 only for some sets of values of b and 
5. Although we do not present here all the energy eigenvalues considered by Adamowski, 
our calculation reveals certain interesting features of this approximation method. 

The present calculations show that the binding spectra of the Hellmann potential possess 
the following properties. 

(i) For low and strong coupling of b in Yukawa part, the energy eigenvalues obtained 
from the perturbation method are in good agreement with the variational results for low 
values of the screening parameter 6. Obviously, when 6 is small the Coulomb field character 
prevails and the perturbation method has been adjusted, of course, to that. However, the 
results become gradually worse as the screening parameter 5 becomes large. Appreciable 
discrepancy of our results from the variational calculations occurs almost for most states if 
5 > 0.2. We suspect that this happen because the perturbative parts of potential becomes 
so shallow and its minimum shifts appreciably from the minimum of the true potential. 
For certain values of 6, the perturbation potentials for some states becomes so shallow that 
the expansion for the energy seies becomes divergent in the sense that higher perturbation 
terms in Eq.(35) dominate over the unperturbed term in Eq.(18) and consequently one gets 
anomalous results. 

(ii) For almost most strongly attractive (negative) b but small S, it is possible to determine 
the binding energy eigenvalues for Is up to 4/ states. 

(iii) For a given n, the results of the energy eigenvalues, Ene increase with increasing i if 
the Yukawa potential is attractive (6 < 0), and Eni decrease with increasing £ if the Yukawa 
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potential is repulsive (6 > 0); i.e., for £ > /, Ene — -B^/ > or Ene — -B^/ < if 6 < or 
b > 0, respectively. This is found to be consistent with the level ordering theorem of Grosse 
and Martin [27]. The s levels are mostly split off from the hydrogenlike levels Ej^ [Eq.(18)] 
(downwards for 6 < and upwards for b > 0). The energy eigenvalues E^e approach E^ if 
n and £ increase. The shift of the energy levels with respect to E^ is due to an influence of 
the finite range Yukawa potential in Eq.(2). 

(iv) For the attractive Yukawa potential there exist some values of b and 6 for which the 
energy levels with larger n and smaller i become lower than the levels with smaller n and 
larger £, i.e., E^e < if n > n' > 3 and i < i' . 

(v) For the repulsive Yukawa potential there are some values of the strength b and the 
screening parameter 5 for which the energy eigenvalues for larger n and i become lower than 
those for smaller n and £, i.e., E^e < if n > n' > 2 and £ > £' . 

On the other hand, in Tables 6 and 7, we present the ground and excited energies (in 
units such that /2m — 1) for several states £ — 0,1,2 calculated for the potential in Eq.(2) 
and compare them with other works [28]. One finds that our results are remarkably good 
for the £ — ground state. As £ increases for a given 6, the error increases. However as b 
increases, the relative error decreases rapidly. This is particularly useful when b is large. For 
large b, the calculated energies compare very well with high precision numerical calculation 
presented by other works [11,12,28]. 

V. CONCLUSIONS 

The detailed analysis of the results in terms of various domains of parameters b and S 
of the Hellmann potential reveals a few important facts concerning the application of the 
perturbed formalism. In the present study the discrete energy eigenvalues for the Hellmann 
potential have been calculated as functions of the strength b and the screening parameter 
S of the Yukawa potential. For 6 = the spectrum is given by Eq. (18). The energy 
eigenvalues Eni for the Hellmann potential are shifted upwards or downwards with respect 
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to the hydrogenlike levels, Eq. (18), if 6 > or 6 < 0, respectively. The absolute values of 
deviations of E^i from decrease with increasing quantum number i, which results in the 
corresponding ordering of the energy levels for a given n. This is due to the influence of the 
finite range Yukawa potential (YP), which decreases with increasing £. 

The properties of the energy spectrum for the Hellmann potential obtained in the present 
work have many analogies in atomic, solid-state, and quark physics. Property (i), giving the 
order of the energy levels Ene with the same n, dependent on £ and on the sign of the YP, has 
an application to such systems as an exciton and a bound polaron in polar semiconductors 
and ionic crystals. These systems consist of two oppositely charged particles interacting with 
themselves through a polarizable medium. The energy levels of both the systems exhibit 
this case, although the total effective potential is even more complicated than the Hellmann 
potential, being a linear combination of the Coulomb potential and an additional potential, 
which is a sum of two Yukawa potentials with different strengths and screening parameters, 
and an exponential potential. However, the net contribution of the additional potential is 
negative for the exciton and positive for the bound polaron and one of the Yukawa potentials 
dominates at small distances [29,30,31]. 

Another system, having the energy levels ordered similarly to those for the Hellmann 
potential with 6 > 0, is the quarkonium (the bound state of heavy quark-antiquark pair). 
For each n the energy levels of the quarkonium system increase with increasing i. This can 
be explained in the frames of the simple model, assuming that the quark interact via the 
potential being the superposition of the attractive Coulomb potential and the positive linear 
potential (Cornell potential) [32]. 
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TABLES 

TABLE I. Energy, —E^e, of states Is — 4/ for Hellmann potential as a function of screening 

parameter 5. It is expressed in units of iJ,a^/2h^ , length in units of oq = Ti^ / jia, and is the 
reduced mass. 



State\(5 0.001 


0.005 


0.01 


0.05 


0.001 


0.005 


0.01 


0.05 


0.1 




T 1 -1 

6 = +1 










b 


= —1 


















Is 


0.250999 


0.254963 


0.259852 






2. 


,24900 


2. 


,24501 


2. 


,24005 


2. 


,20122 






2s 


0.063494 


0.0673535 


0.0719279 






0. 


,561502 


0, 


,55755 


0, 


,552697 


0. 


,517151 






2p 


0.063495 


0.0673775 


0.0720203 






0, 


,561502 


0, 


,557541 


0, 


,552664 


0. 


,516412 






3s 


0.0287644 


0.0324574 


0.0365575 






0, 


,249004 


0, 


,245111 


0, 


,240435 


0. 


,209717 






3p 


0.0287654 


0.0324803 


0.0366436 






0, 


,249004 


0, 


,245103 


0, 


,240404 


0. 


,209055 






Sd 


0.0287674 


0.0325260 


0.0368131 






0, 


,249003 


0, 


,245086 


0, 


,240341 


0. 


,207756 






4s 


0.0166014 


0.0200766 


0.0235509 






0, 


,139633 


0, 


,135819 


0, 


,131381 


0. 


,106316 






4p 


0.0166024 


0.0200984 


0.0236407 






0, 


,139633 


0, 


,135811 


0, 


,131351 


0. 


,105734 
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0.0166043 


0.0201416 


0.0238136 






0, 


,139632 


0, 


,135796 


0, 


,131290 


0. 


,104603 






4/ 


0.0166072 
b = —2 


0.0202058 
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0, 
b 
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= —4 


0, 
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0, 


,131200 


0. 
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Is 
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3.98007 


3. 
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8, 
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8, 
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8, 
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8. 
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8, 
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2s 
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0.990075 
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0. 
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2, 


,24600 


2, 


,23010 


2, 
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2. 


,05964 


1, 


.88728 


2d 


998002 


0.990062 


0.980248 


0. 


90596 


2. 


,24600 


2. 


23008 


2. 


,21033 


2. 


,05807 


1, 


.88135 


3s 


0.442451 


0.434611 


0.425103 


0. 


,359536 


0, 


,996009 


0, 


,980220 


0, 


,960885 


0. 


,820847 


0, 


,678082 


3p 


0.442451 


0.434599 


0.425055 


0. 


,358495 


0, 


,996008 


0, 


,980207 


0, 


,960820 


0. 


,819381 


0, 


.672826 


3d 


0.442450 


0.434574 


0.424959 


0. 


,356435 


0, 


,996007 


0, 


,980174 


0, 


,960691 


0. 


,816464 


0, 


.662442 


4s 


0.248012 


0.240294 


0.231150 


0. 


,174917 


0, 


,558516 


0, 


,542894 


0, 


,524055 


0. 


,397745 


0, 


,288281 


Ap 


0.248011 


0.240281 


0.231103 


0. 


,173993 


0, 


,558515 


0, 


,542878 


0, 


,523991 


0. 


,396397 


0, 


.284004 


Ad 


0.248010 


0.240257 


0.231011 


0. 


,172169 


0, 


,558514 


0, 


,542845 


0, 


,523865 


0. 


,393717 


0, 


.275519 


4/ 


0.248009 


0.240221 


0.230872 


0. 


,169493 


0, 


,558512 


0, 


,542797 


0, 


,523674 


0. 


,389737 


0, 


.262966 



19 



TABLE II. Energy eigenvalues, —Eni, of states Is — 4/ for Hellmann potential as a function of 



scrcxiiiug 


]")arMm('tcr (). Ener 


gy cigcimihu'h 


1 i\rc gi\'c}ii in 


uuils (^f TmV)1(; 1. 






State\(^ 


0.001 


0.005 


0.01 


0.05 


0.1 


0.2 


6 = -10 














Is 


35.99 


35.95 


35.9001 


35.5031 


35.0124 


34.0489 


2s 


8.99 


8.95012 


8.90050 


8.51527 


8.0482 


7.18656 


2p 


8.99 


8.95010 


8.90042 


8.51025 


8.04037 


7.15688 


3s 


3.99001 


3.95028 


3.90112 


3.52703 


3.10435 


2.39157 


3p 


3.99001 


3.95026 


3.90103 


3.52508 


3.09701 


2.36532 


3d 


3.99001 


3.95022 


3.90087 


3.52120 


3.08240 


2.31330 


4s 


2.24002 


2.2005 


2.15197 


1.79676 


1.42667 


0.88467 


4p 


2.24002 


2.20048 


2.15189 


1.79490 


1.41993 


0.86351 


Ad 


2.24002 


2.20043 


2.15173 


1.79119 


1.40652 


0.82123 


4/ 


2.24001 


2.20037 


2.15148 


1.78565 


1.38656 


0.757901 
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TABLE III. Energy eigenvalues,— of states Is — Sd as function of the strength b of the 



Hc'Umaim 


polcnlial for 6 


= (J.Ol. EiK;rgy 


eigenvalueti 


Mix; given in uuil; 


s of Tal)k} 1. 




6\State 


Is 


2s 


2p 


3s 




3d 


1 


0.259852 


0.071928 


0.0720203 


0.036557 


0.0366436 


0.0368131 


0.5 


0.567450 


0.145431 


0.145463 


0.067079 


0.0671090 


0.0671683 


0.2 


0.811983 


0.204435 


0.204446 


0.091858 


0.0918682 


0.0918884 





1 


0.25 


0.25 


0.111 


0.111 


0.111 


-0.2 


1.20801 


0.300553 


0.300545 


0.132562 


0.132553 


0.132537 


-0.5 


1.55753 


0.385743 


0.385723 


0.168871 


0.168852 


0.168815 


-1 


2.24005 


0.552697 


0.552664 


0.240435 


0.240404 


0.240341 


-2 


3.98007 


0.980297 


0.980248 


0.425103 


0.425055 


0.424959 


-5 


12.2 


3.01293 


3.01286 


1.31206 


1.31199 


1.31185 


-10 


35.9 


8.9005 


8.90042 


3.90112 


3.90103 


3.90087 


-20 


120.8 


30.051 


30.0505 


13.2460 


13.2456 


13.2454 
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TABLE IV. Energy eigenvalues, —Ent, of states 5s — 7i as function of the strength b of the 



HcUmaiiu 


polcnlial for 6 


= (J.Ol. EiK;rgy 


eigenvalueti arc 


given ill uuil;- 


i of Tal)lc} 1. 




State\6 


+1 


-1 


-2 


-4 


-8 


-10 


5s 


0.0165535 


0.081148 


0.141758 


0.322393 


0.922921 


1.34306 


5p 


0.0166884 


0.0811191 


0.141713 


0.322331 


0.922844 


1.34298 


5d 


0.0169435 


0.0810614 


0.141624 


0.322208 


0.922692 


1.34282 


5/ 


0.017289 


0.0809751 


0.141490 


0.322023 


0.922463 


1.34258 


55 
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0.141312 


0.321776 
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1.34225 


6s 


0.0100105 


0.0540995 


0.093579 


0.213386 


0.618604 


0.904359 


6p 


0.0103092 


0.054072 


0.0935360 
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0.618530 


0.904281 


6d 


0.0108763 


0.054017 


0.0934503 


0.213206 


0.618380 


0.904123 


6/ 


0.011651 


0.0539351 


0.0933218 


0.213027 


0.618155 


0.903887 


65 


0.0125417 


0.0538266 


0.0931512 


0.212787 


0.617856 


0.903572 


6^ 


0.013427 


0.0536924 


0.0929390 


0.212489 


0.617482 


0.903178 


7s 


0.0028019 


0.0380167 


0.0648977 


0.148193 


0.435796 


0.640565 


7p 


0.0020524 


0.0379905 


0.0648565 


0.148135 


0.435723 


0.640488 


7d 


0.0006123 


0.0379383 


0.0647744 


0.148018 


0.435576 


0.640333 


7/ 


0.0014009 


0.0378606 


0.0646516 


0.147845 


0.435356 


0.640101 


79 


0.0038107 


0.0377579 


0.0644887 


0.147613 


0.435064 


0.639792 


7h 


0.0063819 


0.0376312 


0.0642864 


0.147325 


0.434698 


0.639405 


7i 


0.0088204 


0.0374816 


0.0640455 


0.146980 


0.434260 


0.638942 
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TABLE V. Energy ei{ 


^envalues, —Eni 


, of states Is 


- 6h for 6 = 0.1. 


Energy eigenvalues are 


glA't^ii ill 


uuilH of TaV)l(; 1. 












State\6 


-5 


-8 


-10 


-20 


-30 


-50 


Is 


11.7605 


24.2118 


35.0124 


119.014 


253.014 


671.014 


2s 


2.60282 


5.49595 


8.0482 


28.3034 


61.0555 


164.057 


2p 


2.59637 


5.48852 


8.04037 


28.2947 


61.0463 


164.048 


3s 


0.946371 


2.0766 


3.10435 


11.562 


25.5672 


70.2385 


3p 


0.94056 


2.06971 


3.09701 


11.5536 


25.5583 


70.2292 


3d 


0.929052 


2.05602 


3.0824 


11.5367 


25.5406 


70.2106 


4s 


0.405183 


0.92834 


1.42667 


5.76553 


13.2138 


37.4734 


4p 


0.400245 


0.922129 


1.41993 


5.75749 


13.2052 


37.4643 


4d 


0.390444 


0.909778 


1.40652 


5.74145 


13.188 


37.4461 


4/ 


0.37593 


0.89143 


1.38656 


5.71747 


13.1624 


37.4189 


5s 


0.185308 


0.44099 


0.700006 


3.14653 


7.56577 


22.3833 


5p 


0.181729 


0.43572 


0.694047 


3.13896 


7.55755 


22.3745 


5d 


0.174538 


0.425204 


0.682168 


3.12385 


7.54112 


22.3568 


5/ 


0.16367 


0.40949 


0.664442 


3.10128 


7.51654 


22.3304 


5g 


0.149027 


0.38865 


0.640982 


3.07135 


7.48388 


22.2952 


6s 


0.099768 


0.213846 


0.348659 


1.78538 


4.56702 


14.2628 


6p 


0.097267 


0.209807 


0.343761 


1.77836 


4.5592 


14.2543 


6d 


0.092169 


0.201681 


0.333949 


1.76435 


4.54358 


14.2372 


6/ 


0.084289 


0.189376 


0.319193 


1.74341 


4.5202 


14.2116 


65 


0.073342 


0.172749 


0.299444 


1.71562 


4.48915 


14.17775 


6h 


0.058954 


0.151614 


0.274642 


1.68109 


4.45051 


14.1351 
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TABLE VI. Some bound state energies, —Ene, for neutral atoms. 



b = -Ze^ State Present Calculated [28] Exact [28] 



-16 



Is 


3.25647 


3.219965 


3.250536 


2p 


0.37172 






3d 


0.003891 






Is 


14.4581 


14.419973 


14.457119 


2s 


2.61756 






2p 


2.58365 


2.433176 


2.583677 


3d 


0.54034 






Is 


60.8590 


60.819285 


60.859039 


2s 


13.0276 


13.027315 


13.03259 


2p 


12.9910 


12.837459 


12.991055 


3s 


4.39692 


4.372037 


4.405697 


3p 


4.36357 


4.348041 


4.388576 


3d 


4.29724 






Is 


139.2590 


139.22008 


139.259362 


2s 


31.4314 


31.431281 


31.43595 


2p 


31.3938 


31.238464 


31.393815 


3s 


11.7014 


11.699808 


11.70926 


3p 


11.6662 


11.665284 


11.683877 


3d 


11.5959 


11.245640 


11.595949 


4s 


5 05037 


5 044185 


5 0590 


4p 


5.01801 


5.013466 


5.05410 


4d 


4.95352 


4.951600 


5.00855 


5s 


2.21653 


2.203309 


2.22372 


5p 


2.18793 


2.177040 


2.241432 


5d 


2.13080 


2.124124 


2.24278 
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TABLE VII. Some bound state energies, —Ent, for neutral atoms. 



b = -Ze^ Is 

-13 39.7088 

-36 316.86 

-79 1544.51 



2s 2p 

8.18765 8.15172 

74.0341 73.9958 

374.5 374.461 



3s 3p 

2.57025 2.53816 

29.3129 29.2763 

158.088 158.05 



4s 4p 

0.82329 0.79663 

13.9316 13.897 

82.6361 82.599 
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